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ABSTRACT 

One of the central features of the last 8 to 10 billion years of cosmic history has been the emergence of a well- 
populated red sequence of non-star-forming galaxies. A number of models of galaxy formation and evolution 
have been devised to attempt to explain this behavior Most current models require feedback from supermassive 
black holes (AGN feedback) to quench star formation in galaxies in the centers of their dark matter halos (central 
galaxies). Such models make the strong prediction that all quenched central galaxies must have a large supermas- 
sive black hole (and, by association, a prominent bulge component). I show using data from the Sloan Digital 
Sky Survey that the observations are consistent with this prediction. Over 99.5% of red sequence galaxies with 
stellar masses in excess of IO^^Mq have a prominent bulge component (as defined by having a Sersic index n 
above 1.5). Those very rare red sequence central galaxies with little or no bulge (n < 1.5) usually have detectable 
star formation or AGN activity; the fraction of truly quenched bulgeless central galaxies is < 0.1% of the total 
red sequence population. I conclude that a bulge, and by implication a supermassive black hole, is an absolute 
requirement for full quenching of star formation in central galaxies. This is in agreement with the most basic 
prediction of the AGN feedback paradigm. 

Subject headings: galaxies: evolution — galaxies: bulges — galaxies: spiral — galaxies: general — galaxies: 
stellar content 



1. INTRODUCTION 

A key feature of the last 8-10 billion years of galaxy evolu- 
tion is the emergence of a well-populated red sequence of non- 
star forming galaxi es (e.g.. lBell et"ani2004l : iFaber et al.ll2007l: 
iBrown et alj 120070 . The present-day red sequence has little 
scatter in color at a given magnitude, and appears to be largely 
composed of bulge-dominated and early-type spir al galaxies 
JVisvanathan & Sandagelll977l : lBower et alj|1992l) . Theoret- 
ical models of galaxy evolution have had difficulty in repro- 
ducing the properties of red sequence galaxies using standard 
recipes for the cooling of g as, star formation and stellar feed- 
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back (see. e.g..lBower et al.ll2006;.Croton et al...2006; Cattaneo 



et al. l2007l for discussions of this issue). In such 'standard' 
models, star formation is insufficiently quenched, leading to a 
dramatic overabundance of massive blue star-forming galaxies. 

Current models have attempted to remedy this shortcom- 
ing by quenching star formation in galaxies through two main 
classes of mechanism. The first mechanism, which affects only 
satellite galaxies (i.e., the non-central galaxies in groups and 
clusters), is a shutoff of gas cooling in galaxies once they fall 
into a larger halo. Such a recipe has been in at least semi- 
analytic models for some time (see, e.g.. lCole et alj|200 0) and 
produces a population of relatively low-mass highly-clustered 
red sequence galaxies. These prescriptions are b eing currently 
tested using the SDSS and other surveys (e.g., Baldrv et al.l 
l2006llWeinmann et al.ll200dlHaines et al.ll2007ah . and I will 
not touch further on this issue in this paper. 

Another mechanism (or mechanisms) has been required by 
the models to shut off star formation in central galaxies (i.e., 
the galaxies lying in the center of their dark matter halos, and 
presumably at the focus of any inflow of gas). An important 
candidate (I will touch upon others later) is feedback from ac- 
cretion of matter onto supermassive black holes, which either 
disrupts g as cooling in t he galaxy ha lo (e.g., 'r adio mode AGN 
feedback'. ICroton et al.ll2006., .Fabianet al.ll200 6) or drives gas 
out of the central galaxy (e.g., 'quasar mode AGN feedback'; 



iKauffmann & Haehneltl 120001; iHopkins et al.ll2006l Tremonti 
et al. 1200711 'I In this picture, given the existence of a tight cor- 
relation between black hole mass and bulge mass (e.g. . Magor- 
rian et al. 1998; Hiiring & Rix 2004), it is natural to expect that 
quenching and the existence of a large bulge would be tightly 
linked. 

Indeed, such a corr elation is well-docum ented, at least in a 
broad sense (see, e.g., Strateva et al. 2001 or Fig. 7 of Blanton 
et al. 2003b). Yet, if AGN feedback is the only way for central 
galaxies to quench their star formation, such a paradigm makes 
a strong prediction: bulgeless quenched central galaxies cannot 
exist. In this sense, the small but non-negligible population of 
bulgeless (low Sersic index) galaxies with red colors in Fig. 7 of 
Blanton et al. (2003b) is of key importance. Are these galaxies 
all satellite galaxies and/or dust-obscured edge-on galaxies? Or, 
is there indeed a significant population of red, central bulgeless 
disk galaxies? In the latter case, one would be driven to prefer, 
at least in certain circumstances, other possible mechanisms for 
quenchin g star formation (see. e.g...Naab et al. 2007. Dekel & 
Birnboiml2007L iKhochfar & OstrikeHl2007l iD ekel & Birnboinj 
l2006llBirnboim et al.ll2007l or lGuo & Ohll200 7 on gravitational 
heating, the influence of the development of virial shocks, and 
the heating of large halos with cosmic ray energy). 

The goal of this paper is to attempt to test the key predic- 
tion of AGN feedback — namely, that bulgeless quenched cen- 
tral galaxies do not exist — using the Sloan Digital Sky Sur- 
vey's Data Release 2 ( Abazajian et al. 2004). In ^ I describe 
the data and derived parameters that I use in this investiga- 
tion. In ^ I describe the results, showing the properties of 
the galaxy population as a function of galaxy morphology for 
central galaxies. In sections |4] and |5] I discuss the results and 
present my conclusion. I adopt//o = 70kms~' Mpc~', 51a = 0.7 

'it is expected that such mechanisms act together; e.g., stellar or quasar- 
powered winds may expel the gas from the galaxy originally, while low-level 
AGN activity may suppress future gas cooling. Most models include both types 
of feedback. 
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Fig . 1 . — The color-mass distribution for central galaxies with h/a > 0.5 and 0.02 < z < 0.06. Galaxies are separated into those lacking significant bulges (Sersic 
index n < 1.5; left), those with intermediate bulge-to-disk ratio (1.5 < n < 3), and those with significant bulges (n > 3; right). Colors denote the emission-line 
classification classes from Brinc hmann et alj J2b04) . Red denotes galaxies unclassified because of a lack of line emission; the bulk of these galaxies are not forming 
stars. Hues of blue show galaxies with line emission characteristic of being powered by star formation; green shades show star formation/ AGN composites or AGN. 
There are very few m < 1.5 central red sequence galaxies, and almost all of these have an AGN. It would appear that in order to quench star formation completely, 
central galaxies must have a significant bulge and/or supermassive black hole. 



and r^n, = 0.3 in what follows. 

2. THE DATA 

I use publicly-available catalogs derived from the Sloan Dig- 
ital S ky Survey's (SDSS) Data Release Two (Abazaiian et al. 
l2004(n . The SDSS is an imaging and spectroscopic survey that 
has so far mapped ~ 1 /4 of the sky. Imaging data are produced 
si multaneously in five photometric bands, namely u, g, r, i, and 
zjFukugita et al.l ll996t lOunn et al.l [T998';'Hogg et all 1200 lb 
iGunn et al. 2006) . The data are processed through pipe lines to 
measure photometric and astrometric properti es (Lupto net al.l 
T999l:IStoughton et al.ll200llSm ith et al ." 20021: .Pier et al.ll2003b 
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Ivezic et ai]|2004l : lTucker et aL. 2006) a nd to se l ect targets for 



spect roscopic follow-up (IB Ian ton et al.l 12003d : IStrauss et al.l 
l2002h . DR2 includes spectroscopy over an area of ~ 2600 
square degrees, and imaging for a larger area. 

For this paper, I use the sample of galaxies in the New York 
University Value Added Galaxy Catalo g (NYU VAGC : Blan- 
ton et al. 1200 5). For galaxies with spectroscopic redshifts, I use 
absolut e magnitudes A:-correc ted to the rest-frame g and r pass- 
bands dBlanton et al.ll2003a[) and corrected for G alactic fore- 
ground extinction following 'Schleg el et al.l (Il998l) . I adopt the 
r absolute Petrosian magnitude for the galaxy absolute mag- 
nitude (random and systematic uncertainties < 0.15 mag), and 
the {g- r) model color for a higher S/N estimate of galaxy color 
(uncertainties < 0.05 mag). Stellar masses were estimat ed un- 
der the assumption of a universally-applicable Chabrieij (l2003h 
stellar IMF using the following color-stellar M/L relation: 

logioM/U = -0.406+ 1.097(^-r), (1) 

wher e an offset of -0 .1 dex has been applied to the relation 
f romlBell et al.l (l2003l) to convert from a diet Salpeter IMF to 
a lChabrien (l2003h IMF. These stellar masses have random un- 



certainties of r^ 0.1 dex: systematic uncertainties from recent 
bursts of star formation and other sources may exceed this esti- 
mate in cases where the bursts contribute a significant fraction 
of the galaxy's fight (Bell & de Jong 2001) . 

^DR2 is used in preference to larger, later SDSS data releases because it is 
the only data release to have publicly-available galaxy density estimates, group 
catalogs, Sersic fits and star formation rate estimates. 



In order to gain a more complete understanding of star for- 
mation and AGN activity in the sample galaxies, star formation 
and AGN classificati ons and estimates of to tal star formation 
rate were taken from Brinchmann et al. (2004), using emission 
line measurements described in Tremonti et al. (2004). Galax- 
ies are classified as star forming, AGN, composites, or are left 
unclassified (typically because the galaxies lack line emission 
in their SDSS spectra). 

Disk galaxies were select ed by selecting galaxies with nearly 
an exponential light profile. iBlanton et al.l l2003b) fit th e light 
profile of galaxies in the SDSS with a seeing-convolved feersicl 
(119681) profile E(r) ex exp[-K.{j-^ '/"- 1)]. The Sersic index n de- 
scribes the shape of the light profile, where n=\ corresponds to 
an exponential light profile and n = 4 corresponds to a r'/"* law 
profile characteristic of massive early-type galaxies. In what 
follows, I apply a cut to the Sersic indices in r-band n < 1.5 to 
select disk galaxies with little or no contribution from a bulge to 
the light profile of the galaxy; I apply a further cut that the half- 
light radius should be larger than 0.4" to discard all galaxies 
with litde/no structural information in the NYU VAGCQ. 

Finally, in order to test the model predictions for the proper- 
ties of central and satellite galaxies separately, I use a volume- 
limited catalog of galaxy groups complete to log, „ Mgrpu p iMp^ > 
11.7 for group redshifts z < 0.06 from lYang et alj ( l2005h . 
Groups are identified using an iterative method which adopts a 
trial mass and identifies galaxies within a group mass-dependent 
transverse radius and velocity difference. Then based on the to- 
tal r-band luminosity of all candidate group members the group 
mass is adjusted and the process repeated until one converges 
on a final group mass. The group mass is assigned based on the 
total r-band group luminosity of all members; in a cosmolog- 
ical simulation with the same volume as the observed sample, 
the most luminous group is assigned the largest halo mass, the 
second most luminous group the second largest mass, etc. Such 
groups are ^ 90% complete and contaminated at the ^ 20% 

^The Sersic fits are used in preference to the more frequently used concen- 
tration parameter c,- because the Sersic fits are convolved with the point spread 
function (i.e., the Sersic fits should be less dependent on seeing than the seeing- 
dependent concentrations). 
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Fig. 2. — SDSS cutouts of galaxies in tlie sample with M, > W'^Mq. Top: Blue n < 1.5 central galaxies. Middle: Red n < 1.5 central galaxies. Bottom: Red 
1.5 < n < 2.5 central galaxies. The rightmost galaxy in the middle row is the only red n < 1.5 central galaxy lacking Hne emission. Cutouts are roughly 40" on 
a side. There is a tendency for the red n < 1.5 central galaxies to appear rather more concentrated than their blue counterparts, lending weight to the notion that a 
significant bulge component is a requirement for quenching of star formation in central galaxies. 



level. In this group catalog, the majority of galaxies reside 
in single galaxy 'groups'; i.e., they are the central galaxies in 
their own halos, with no satellites above the luminosity limit of 
the sample. For the purposes of this paper, I am almost com- 
pletely unaffected by group mass uncertainties, as the primary 
purpose in using the catalogs is to separate the galaxy popu- 
lation into central galaxies and satellites; obviously, for single 
galaxy groups the identification of the central galaxy is trivial. 
For all multi-galaxy groups, Yang et al. identify the brightest 
galaxy in the group as the central galaxy. Comparison with 
mock group catalogs derived using this method demonstrates 
that > 97.5% of central galaxies are in deed the brightest galaxy 
in their group JWeinmann et al.ll2006l) . Satellites are defined to 
be all other galaxies in groups. 

3. RESULTS 

The goal of this paper is to test the generic prediction of AGN 
feedback that there should be no bulgeless central galaxies with 
quenched star formation. In order to test this prediction, I select 
galaxies with 0.02 < z < 0.06 from the above sample. I further 
select galaxies with < 60° inclination (b/a > 0.5). This is a 
selection cut of some importance, as it minimizes the effects of 
dust extinction on galaxy colors and emission line-derived esti- 
mates of star formation rates; low-inclination galaxies with red 
colors and/or low star formation rates become then a reasonably 
clean probe of quenchingj. 

'^A further advantage of the selection of low-inclination systems is that the 
circular apertures used for the SDSS model and Sersic fits are more appropriate 
for these more cu'cular systems than they are for higher inclination systems. 



The key result of this study is shown in Fig. [T] The left-hand 
panel shows the properties of central galaxies without signifi- 
cant bulges, selected by requiring that the Sersic index n < 1 .3^ 
The right-hand panel shows the properties of central galaxies 
with large bulges, selected using n > 3. It is to be noted that 
the n > 3 bin contains a number of galaxies with prominent 
(even dominant) disks; the « > 3 is reflecting primarily the exis- 
tence of a large bulge. The middle panel shows the intermediate 
1.5 < n < 3 bin. One can see that there are very few bulgeless 
galaxies on the red sequence; adopting a definition of red se- 
quence of (g - r) > O.57 + O.O5751ogio(M,/lO**M0) (i.e., aline 
0.05 mag bluewards of the red sequence locus described by the 
thick line) and restricting our attention only to galaxies with 
logioM»/M0 > 10, only 0.5% (14/2744) of red sequence cen- 
tral galaxies have « < 1.5 (14/725, or 1.9%, of n < 1.5 galaxies 
are on the red sequence jj. 

^The referee pointed out that there are a number of examples of low-mass 
red sequence early-type (E/SO) galaxies with n < 1.5, e.g., in the Virgo cluster. 
Indeed, in the current sample, the vast majority of the red, non-star-forming 
n < 1.5 galaxies are low-mass (M» < W^'^Mq) satellites, primarily of groups 
with group masses in excess of IO'^^Mq. There are no « < 1.5 non-star-forming 
central galaxies with M, < lO'^M© (Fig.[T] left panel), lending weight to the 
notion that while non-star-forming n < 1.5 E/SO galaxies exist, they exist al- 
most exclusively as satellites of reasonably massive groups, and are the product 
of (primarily hydrodynamical) stripping of late-type disk-dominated galaxies. 

^There are 16 red sequence central galaxies with values of h/a > 
0.5 and n < 1.5; two of these (SDSS J142922.73-004939.1 and SDSS 
J084958. 78+381203. 2) are clearly edge-on and should have not been included 
in the sample (their catalog values of b/a are simply incorrect) and they were 
excluded from the sample. 
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Fig. |2] shows example gri SDSS cutouts for subsamples of 
central galaxies with M* > IQ^'^Mq. The top panels show five 
of the 711 blue central n < 1 .5 galaxies; one can clearly see that 
the typical galaxy in this subsample is a late-type disk galaxy 
with a small or nonexistent bulge. The middle panels show five 
of the 14 red central n < 1.5 galaxies. One can see in some 
cases a prominent bar and possible bulge. The lower panels 
show five of the 163 red central galaxies with 1.5 < n < 2.5. 
There is a tendency towards better defined bulges than for the 
n < 1 .5 red central galaxy subsample; however, it is clear that at 
least some fraction of the low-n red subsample have structures 
similar to the higher n subsample, lending weight to the notion 
that some of the 14 n < 1 .5 red central galaxies have incorrect n 
estimates — 0.5% is an upper limit to the fraction of bulgeless 
red sequence galaxies. Taking Figs.[T]and|2]together, it is clear 
that there are very few, if any, bona fide bulgeless red central 
galaxies — in agreement with the most basic expectation of the 
AGN feedback paradigm. 

One may have the concern that perhaps this result was an 
artifact of the way in which central and satellite galaxies are 
chosen by the group finding algorithm. One can repeat the 
analysis using cylindrical overdensities from the NYU VAGC 
(i.e., overdensities measured in cylinders of 1 .4 Mpc radius an d 
1600kms"' extent along the line of sight; jBlanton et all2003bl) . 
which are simpler and more conventional measures of environ- 
ment. Isolated galaxies (defined as having overdensities (5 < 0) 
are likely to reside as centrals in their own halos. When re- 
peating the analysis using isolated galaxies as a (conceptually 
less optimal) proxy for central galaxies, one finds that the frac- 
tion of M* > 1O'"M0 red sequence galaxies in isolated envi- 
ronments with n < 1.5 is 1.0% (18/1842). So one finds, us- 
ing a completely different technique for identifying likely cen- 
tral galaxies, the same result: in central galaxies, quenching is 
empirically correlated with the existence of a prominent bulge 
component. 

4. DISCUSSION 

It is worth asking if this empirical association between hav- 
ing a bulge and the ability to quench is a truly bulge-specific 
association, or if it reflects a third, unexplored parameter. An 
obvious candidate parameter is surface densitv. In Kauffmann 
et al.l2003) and .Bell & de Jong. (2000.) it was found that the 
stellar surface density of a galaxy was correlated strongly with 
star formation history, inasmuch as galaxies with high surface 
densities appeared to form the bulk of their stars at earlier 
epochs than galaxies with lower surface densities. In Fig. [3] 
I plot optical g-r rest-frame color as a function of the stel- 
lar surface density /i* within the half light radius for bulgeless 
galaxies (« < 1.5; black) and galaxies with large bulges (n > 3; 
gray). While g-r is indeed a function of surface density for the 
n < 1.5 galaxies, it is clear that at a given surface density (e.g., 
in the range around ^ lO^Mokpc"^) that galaxies with strong 
bulges are primarily quenched, whereas bulgeless galaxies at 
that same density still actively form stars. Such a conclusion 
was also reached by Kauffmann et al. (2006, see their Fig. 14), 
who explored this star formation rate per unit stellar mass as a 
function of surface density for a samples of galaxy split by con- 
centration. Thus, one can rule out that the correlation seen in 
Fig. [T] is driven by underlying correlations with stellar surface 
density. 

Fig.|4]shows the distribution of the Sersic indices of the sam- 
ple of central galaxies. One can see that the fraction of the 
sample that is red is a strong function of the Sersic index, with 
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Fig. 3. — The color-surface density distribution for central galaxies with 
h/a > 0.5 and 0.02 < z < 0.06. Galaxies with large bulges (n > 3) are shown 
in gray, and bulgeless galaxies (n < 1.5) are shown in black. Quenching is 
not driven by surface density; at all surface densities, galaxies with prominent 
bulges n > 3 have much redder colors than n < 1.5 systems. 
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Fig. 4. — The upper panel shows the histogram of Sersic values for all central 
galaxies (gray), red central galaxies (dark gray), and galaxies with no detected 
emission lines (black). In the lower panel, the fraction of central galaxies that 
is red (dark gray) and that have no emission lines (black) is shown as a function 
of Sersic index. 



a small red fraction for n < 2 smoothly increasing to n ~ 4.1 

'The red fraction does not asymptotically approach unity for the high n 
systems because of the combined influence of our somewhat restrictive red se- 
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Interestingly, only a small fraction (^ 1/4) of red galaxies have 
no emission lines at all (the fraction of galaxies which are com- 
pletely quenched, at least to the detection limits characteristic 
of the SDSS; see Haines et al. 2007b for the same result using 
GALEX data in conjunction with the SDSS). Put differently, 
complete quenching of gas infall into in central galaxies is a 
relatively uncommon occurrence; most red sequence centrals 
have some modest AGN activity or low-level star formation. In 
particular, only 1 out of the 14 red central « < 1 .5 galaxies with 
stellar masses above lO^'^M© lacks line emission Fig. [1] i.e., 
complete quenching of bulgeless galaxies is a very rare occur- 
rencq^. 

Finally, it is obvious that while bulges appear to be a neces- 
sary condition for the suppression of star formation, the simple 
existence of a bulge is not a sufficient condition to suppress star 
formation. There are a number of suggestions as to what this 
'second parameter' cou ld be: pseudobulges vs. classical bulges 
JDrorv & Fisheij|2007h . th e transition to a hot viri alized halo 
at h alo masses ^ IO'^M^t^ JDeke l & Birnboimll2006,; Cattaneo 
et al. l2006l) , or the rat e at whichc osmological infal l can deposit 
energ y into the halo dNaab et all 120071: iKhochfar & Ostriked 
l2007h . In this context, it is noteworthy that for group masses 
> IO'-'Mq 23% of the n > 3 central galaxies are blu^l (as op- 
posed to 35% for halos with masses below lO'^M©). Thus al- 
though models which assume quenching is co mplete above a 
give n mass cut (e.g., lO'^M© in the work of Cattane o et al.l 
l2006i) offer a qualitatively interesting and computationally con- 
venient approximation for the evolution of the galaxy popula- 
tion, it is clear that such an approximation is not correct in detail 
— even at group masses > 1O'^M0 quenching is incomplete. 

5. CONCLUDING REMARKS 

The main message of this paper is simple. Galaxy evolu- 
tion models with AGN feedback predict that quenching of star 
formation in central galaxies should occur only in those cen- 
tral galaxies with prominent bulges (therefore, large supermas- 
sive black holes). I have shown that this prediction was essen- 
tially correct, using data from the SDSS DR2. I choose to an- 
alyze only relatively low-inclination galaxies with b/a > 0.5; 
while the general conclusions of this paper would hold also 
if higher inclination galaxies were included, the scattering of 
n < 1.5 galaxies onto the red sequence by dust extinction con- 
siderably dilutes the significance and impact of the results. I 
found that at least 99.5% of red sequence galaxies with stel- 
lar masses > lO'^M© have prominent bulges (Sersic indices 
n > 1 .5). Furthermore, most the 0.5% of bulgeless red sequence 
central galaxies had detectable line emission; there is only 1 
n < 1.5 central galaxy without detectable star formation out of 
the 6036 central galaxies in this mass-limited sample. This al- 
most perfect empirical association between having a prominent 
bulge and quenching is in excellent agreement with the basic 
expectation of AGN feedback. 

I thank the referee for their constructive comments and ques- 

quence cut and the scatter of a small fraction of high n systems towards the 
blue, as seen in Fig.[T] 

*The right-most panel of the middle row of Fig. |2] shows this galaxy; this 
galaxy appears to be significantly more compact and concentrated than its n < 
1.5 would imply. 

'Most of these blue n > 3 central galaxies are spiral galaxies with large 
bulges (70%), 20% appear to be relatively undisturbed early-type (E/SO) galax- 
ies, and 10% of them are highly disturbed, with asymmetries or tidal tails in- 
dicative of a previous interaction or merger. 
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